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ABSTRACT 


The vertical structure of the mean longshore current is examined on three strong 
current days during the DUCK 94 experiment and it is found well described by a logarithmic 
profile (mean correlation coefficient for all 22 profiles, 0.98). This hypothesis works better 
in the trough where turbulent bottom boundary layer processes are predominant than over 
the bar, where breaking-wave induced turbulence generated at the surface modifies the 
profile. 

A relationship between near bottom vertical velocity profiles and bottom roughness 
was found. The bed shear stress coefficient varied by an order of magnitude across the surf 
zone (0.001-0.05). For the three days considered, it is concluded that the bed shear stress 
coefficient increased with increasing bottom roughness, and therefore is an important 
parameter to characterize the bottom boundary layer. 

The influence of the wind is parameterized by an exponential approximation that 
works well in 50% of the analyzed profiles. This parameterization assumes that the residual 
Data-Normal distribution of data is only due to the wind and does not consider the 
alongshore component of mass transport velocity, which is assumed small. 

These results should be treated carefully since there were obtained for only three days 
over a barred beach with strong longshore currents, and cannot be generalized without 


further studies. 
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I. INTRODUCTION 


Knowledge of the bottom and surface boundary layers is fundamental to 
understanding nearshore hydrodynamic and sediment processes. For steady flow, such as 
in a river, the bottom boundary layer is well described by a logarithmic profile. The 
superposition of waves on the mean current produces enhanced bottom friction (e.g., Grant 
and Madsen (1979), Christofferson and Jonsson (1985), Myrhaug and Slaattelid (1989) and 
Sleath (1990 ) ). As a result, the gradient of mean current near the bed is decreased and a 
more linear profile can be expected. The influence by waves is inversely proportional to the 
depth of the water with decreasing importance of wave-bottom boundary with increasing 
depth. 

Only limited wave-driven, longshore current vertical profiles have been measured. 
Visser (1984) measured wave-driven longshore currents in a laboratory experiment using 
dye displacement over the vertical. He found the alongshore spatially averaged currents 
showed more linear than logarithmic profiles. However, Visser (1986), repeated the 
experiment in the same basin, measured the vertical profiles using micro-propeller and laser 
Doppler anemometers and found profiles approached more a logarithmic form. In a similar 
laboratory wave-driven longshore current experiment using laser nemonees Simons ef 
al.(1992) found the vertical velocity profiles tend to show logarithmic behavior. Church ef 
al.(1993) used field data from SUPERDUCK experiment (1986) to study the problem. Three 


electromagnetic current meters mounted on a mobile sled were used in a vertical stack 0.7, 





1.0 and 1.5 m above the bed. The results were not conclusive in relation to the logarithmic 
profile question since only three current meters were used. It was concluded that more 
current meters located nearer the bed are required to determine the bottom boundary layer 
shape. 

Turbulence induced by breaking waves can modify the vertical profile of longshore 
currents. The downward vertical momentum mixing produced from the breaking wave- 
injected turbulence results in a more uniform velocity profile, increasing the shear stress at 
the bed ead et al., 1986 and Church and Thornton, 1993). 

Wind also modifies the vertical profile of longshore currents. Murray (1975) 
derived the vertical profiles of mean currents as a function of the wind in shallow water by 
balancing the surface stress by wind with the bottom stress, with and without the Coriolis 
force. He compared the theory with observations of current profiles obtained using drogues 
at various elevations well outside from surf zone. Comparing wind and wave forcing of 
longshore currents within the surf zone, Whitford and Thornton (1993), showed the wind 
can be significant, but only at high speed and when the wind has a strong alongshore 
component. 

In the following, mean longshore current profiles obtained over a barred beach are 
examined with the following objectives: 

a. Test the hypothesis that the turbulent bottom boundary layer of the mean longshore 
currents is logarithmic. 


b. Determine how breaking-wave induced turbulence generated at the surface 


modifies the vertical velocity profile. 








c. Investigate how wind modifies the vertical profile of the longshore currents. 
d. Examine the relationship between near bottom vertical velocity profiles and 


bottom roughness, including the influence of ripples and mega-ripples. 














Hl. THEORY 


The vertical profile of the longshore currents is significantly affected by the physical 
boundary layers of the bottom and sea surface. The bottom boundary assumes more 
importance since it determines the general logarithmic profile shape. However, processes 


in the surface layer can modify the profile in the presence of waves and wind. 


A. BOTTOM BOUNDARY LAYER 
Neglecting molecular viscous stresses, the alongshore momentum equation ( y- 


direction) is written 


dpv : Opuv : dpv? : dpwy _ _ op 
or Ox oy dz oy 








(1) 


Let u,= U,t+u ,+a, and w= w ,+W., where the horizontal (u,, i=1,2) and vertical w 
velocities have been expanded into mean, turbulent and wave-induced components. Assume 


straight and parallel contours, i.e., <0 (overbar indicates time averaging), neglect 
M4 ar — 
horizontal turbulent and mean momentum mixing * _ a and neglect = since its 
x x Zz 








contribution is small above the wave bottom boundary layer. Assuming steady state 


conditions and after time averaging 


Opvii _ _ dpw'v 
Ox Oz 





(2) 








which says the cross-shore change in the wave induced alongshore momentum flux 1s 
balanced by vertical changes in alongshore turbulent shear stress. In the shallow water of the 
surf zone, vz can be assumed independent of depth, and its gradient in x is constant in the 


y-direction for straight and parallel contours. Defining 
T= -pw'v' = pv — (3) 
where v, is the turbulent eddy viscosity. Substituting (3) into (2) 


Ty _ Opti 
Oz Ox 


= R (4) 
where # is constant alongshore. Constant shear stress distributions occur in flows driven by 
constant hydrostatic pressure gradients, such as in steady open channel flow, which are well 
described by a logarithmic velocity profile. Therefore it is hypothesized that the vertical 
velocity profile for longshore currents 1s logarithmic. A steady, uniform, turbulent boundary 


layer in the alongshore direction can be described by the logarithmic “law of the wall’: 


y 
* in ( zth 
K Z 


oO 








v(z) = ) (5) 


where z is positively upward from the surface, / is the depth of water, x is the Von 
Karmann constant (0.4), v, is the shear stress velocity and z, is the roughness height. The 
effect of including waves on the bottom stress is to effectively increase the apparent 


roughness height denoted by replacing z, with z, (that is analogous to, but larger thanz, 
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since it includes the additional stress caused by the waves). Following the empirical model 
by Nielsen (1992) for current profiles in the presence of waves, z_ can be related to the 


thickness of wave-dominated layer 6 by 


z =e b6 (6) 


Since 6 has no generalized formulation, a particular formula is suggested. For fairly rough 


conditions (0.06 < m < 0.5 ), constant eddy viscosity, and neglecting the angle between 


wave and current, 6 can be parameterized as 


5 ~ 0.016 Aw 
. . (7) 


K a 








where @ is the radian frequency, r is the hydraulic roughness and A is the orbital amplitude 
of the waves just above the bottom wave boundary layer. 
The shear stress profile is determined by integrating (4) and assuming the stress at 


the surface is equal to zero, to give 
t(z) = -t(-h) 2 8 
yi) = ~t(-h) h (8) 
where / is the total depth. Bottom stress is related to v, through 


t (=~h) = ov? (9) 


Substituting (5) and (8) into (3) we obtain 








v= “hy 21+) - (10) 


which says the eddy viscosity is parabolic with depth. However, since it is based on the “law 
of the wall”, (10) would not be expected to be an adequate formulation near the surface, 


where the wave and wind effects can be significant. 


B. SURFACE BOUNDARY LAYER 

The surface boundary layer is modified by the waves and wind. Wave modifications 
can be due to three contributions: modification of mean currents by measuring currents in 
a Eulerian frame with an undulating boundary; contribution to longshore currents of mass 
transport velocity in the crest-trough region by obliquely incident waves in a Eulerian frame, 
and modification by the wave breaking. 

The most significant modification of the near surface mean current profile is due to 
the effect of hand currents in an undulatory wave boundary in a Eulerian frame. In the 
wave crest/trough region, the averaged measured velocity is decreased because the current 
meter is out of the water part of the time. For instance, the current meter is out of the water 
half the time at the MSL, and the time averaged current is only 50% the expected value of 
the logarithmic profile. To account for this, the surface elevation probability distribution 
function (pdf) is applied to the expected mean current profiles in the absence of waves. The 
percent of time the current meter is in the water is given by 1-P(n), where P(1) is the 
cumulative surface elevation pdf In a Eulerian frame of reference, the modified mean 


current in the crest-trough region is given by 











Vi (2) = [1-P(n V2) (11) 

where V(z) represents the log profile in the absence of waves. For mild wave Heielits and 
in deeper water, the surface elevation pdf is well described by the Gaussian pdf As will be 
seen, the measured pdf's in the surf zone are slightly positively skewed from the Gaussian 
distribution. 

The mass transport velocity can contribute to the longshore current in the upper 
boundary layer when the incident wave angle is relatively large. Assuming irrotational flow 
and evaluating mass flux (transport), AZ by considering separately two regions in a Eulerian 


frame of reference (cf. Philips, 1977) 





M=p (f° ,adz + {judz) (12) 


The first term is zero for irrotational flow. For infinitesimal linear wave theory, wz is not 
defined between zero and the surface n. A Taylor series expansion about z=0 1s used to 


extend defined values of u(0) to the surface, giving 


M = pu(0)n (13) 
correct to O(a/k)”. Applying linear theory 
ee (14) 


where £ is the energy and C the celerity. The mass transport for a single wave is interpreted, 
to this order in an Eulerian reference frame, as due to a uniform velocity confined to the 
crest/trough region. The mass transport velocity is defined 


9 








M 
ea - (15 
ay (15) 


ui 


where A 1s the amplitude. 

For random waves, the wave amplitudes are reasonably described by the Rayleigh 
distribution. The only waves that contribute at any elevation z, will have an amplitude 
A > z. Hence, the ensemble average mass transport velocity profile is obtained by applying 


the wave amplitude pdf 


<U(z)> = [ _UArp(A)dA (16) 


The Rayleigh distribution is given by 





> e = (17) 





H (- y ft 

W 2Z H Tt 2Z 

<U(z)> = ————[—— e ™ 3+ += erfe(——) (18) 
@) 8 tanh th 2 J | 


rims rms 


where erfc(x) is the complementary error function. The alongshore component is defined by 


<V(z)> = <U(z)>sin0 (19) 


where 6 is the mean wave incident angle. 











The surface boundary layer during wave breaking is typified by an intense 
production of turbulence that eventually is dissipated in the shear layer at the lower 
boundary of the surface roller. The remaining energy propagates vertically through the water 
column. The increase of eddy viscosity due to breaking-induced turbulence produces a more 
vertical profile of the current for a given shear stress distribution, compared with profiles 
in the absence of breaking. This results in a stronger velocity gradient near the bed, 
increasing the bed shear stress. 

The influence by wind on the longshore currents can be significant during times of 
strong alongshore wind. Again integrating (4), but including an alongshore wind stress 


OT 
[ ie dz =) - t,(z) = R (n-z) (20) 


with T= p_C,,WWcos®@ (surface) and a = pC,|u|v, (bottom), where W is the wind speed, p, 
is the density of air, @ the wind angle with the longshore current, C,, the drag coefficient, C : 
the bed shear stress coefficient andv, the longshore current speed at the bottom. Equating 


the stress at the surface with the alongshore wind stress component we have 


ov 
ty a Pv (21) 
Solving for v exclusively due to the wind 
v= vec Weosh|”, (22) 
p V, 


I] 





Murray (1975) considered constant v, 


Ss 
vy = ve cosh h (23) 
p Y, 


This theory applies when the wind angle to the coast is less than 45 degrees and does not 
consider the influence of the waves. It is found this treatment overestimates the speed at the 
bottom, which was verified using field data from DUCK94. 

More complex forms can be adopted for eddy viscosity such as the half-parabolic 


distribution (Roelvink and Reniers, 1994). Herev, is defined as 


2 ae 4 Z 
V=—v— Q-— 24 
with 
_ 7 
see ape (25) 
3 p 


As will be seen, an exponential distribution of the wind-driven velocity appears to 


give a reasonable fit to the data 


v= ve uF (26) 


where v, is the velocity at MSL (top current meter) and B a free parameter. The 


corresponding eddy viscosity is 


12 





—_— 


2 
go Pa OD a 
Pp 


8B (27) 


t 


This approach is used to interpret the data from DUCK94. 
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III. DUCK 94 EXPERIMENT 


The measurements described here are part of the comprehensive nearshore DUCK94 


experiment conducted during August and October 1994 at the U.S. Army Corps of 
Engineers Field Research Facility (FRF), Duck, North Carolina. The FRF is located on the 
Outer Banks, a barrier island formation with no major coastal structures that can obstruct 
nearshore flows. The beach is a two-bar system with a dynamic inner bar (30-120m 
offshore) and a secondary bar with lower amplitude (300-400m offshore). The mean 
foreshore slope of the beach is approximately 0.08 (1:12) and the slope offshore of the bars 
is approximately 0.006 (1:170) (Lippmann, 1993). The mean tidal range is 1.0 m. Sediments 
within the surface zone are well sorted with a mean grain size of 0.2 mm. Sediments on the 


foreshore are poorly sorted with larger mean grain size (0.4 mm). 


A specially designed sled was used as a platform to ane instruments (Fig. 1). The 
sled is constructed of a 3 x 4 m six-inch aluminum-pipe frame with two 5m, 8 inch pipe 
runners. This low-profile structure was stabilized by 180 Kg of lead weight plus 
approximately 450 Kg of sand inside the runners. In addition, four fins (45 cm wide) 


extended 60 cm into the sand to insure the sled did not move in the surf zone. 


The vertical profile of longshore current was measured using a vertical stack of 8 
Marsh-McBirney two component electromagnetic current meters with 2.5 cm diameter 


spherical probes mounted on a 2.5 m mast (Fig. 1). The current meter elevations above the 


1S 











bed were 23, 42, 68, 101, 147, 179, 224 and 257 cm. The sled was oriented such that the 
vertical slack of current meters was consistently placed on the “up-current” side erie Sled 
to avoid flow contamination by the sled structure during observations. The current meters 
were pre- and post- calibrated in a tow tank at the Naval Postgraduate School with an 
agreement of 1.9 % in gain. An in situ determination for the offset is used, which was 
obtained by reversing the orientation of the current meters on a very slow longshore current 
day (8 October) by turning the sled around and returning it to the same location (within 1 


m) within one hour. The in situ determined offsets were within 1 cm/s. 


Waves and mean water level were measured using an array of six pressure sensors 
configured in a 3m square with one in the center and another in the middle of a cross-shore 


leg. 


The data were digitally encoded to 14 bit precision on the sled at 36 samples/second 
and transmitted to shore via a fibre-optic cable where signals where monitored and recorded. 
Short cables from the sensors to the data acquisition system on the sled (<7m) resulted in 
exceptionally low noise electromagnetic current meter and pressure sensor signals. An 
armored cable, married to the sled chain tether, provided power and controller signals for 
the instruments via two conductors and returned the digitized signals and video via two 


fibre-optic lines. 


The sled was towed to the farthest offshore location for the first run (approximately 


160m from the shoreline) by the Coastal Research Amphibious Buggy (CRAB). A four- 
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wheel drive forklift pulled the sled shoreward 10 to 30 meters for subsequent runs (each run 
was nominally one hour). Five to eight runs were made across a transect each day. The data 
were acquired during daylight to early night. 

The morphology of the bottom (bathymetry) was measured at various scales from 
the CRAB. Large-scale variations of bathymetry were obtained by using an autotracking 
laser ranging system to measure the CRAB position with a vertical accuracy of less than 3 
cm. Contour plots of the bathymetry for the days selected for analysis show the contours in 
the vicinity of the sled measurements to be quasi straight and parallel (Fig. 2; October 11th 
is typical). Small-scale vertical bottom variations relative to the CRAB, including ripples 
and megaripples, were acquired with a 1 MHz sonic altimeter mounted on the CRAB 70 cm 
from the bed with a vertical accuracy less than 2 cm (Gallengher ef a/., 1995). Bathymetry 
for selected days (Fig. 3) shows a pronounced bar progressively moving offshore and 
significant small-scale morphology in the trough. Areal variations were determined using 
a 500 Khz side-scan sonar also mounted on the CRAB. 

Meteorological information of wind, air temperature, atmospheric pressure and sea 
surface temperature were recorded simultaneously at the end of the FRF pier and onshore 


near the FRF laboratory. 
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IV. DATA RESULTS 


The data analysis presented is for the October phase of the DUCK94 experiment. The 
weather during October is climatologically characterized by three distinct phases (Fig.4): 
weak currents and winds from north (4-9 Oct.), relatively strong currents from north (0.6-1.0 
m/s) caused by a storm with predominant winds and waves from north (10-12 Oct.), and 
variable currents and winds from north/south (13-21 Oct.). The data selected to analyze are 
10-12 October, during the strong longshore currents period. 

All 22 vertical profiles of longshore currents obtained during the three days are used. 
The profiles are based on the measurements by seven current meters over the vertical. The 
current meter the closest to the sea bed was not used because of malfunction. Mean 
longshore velocities are 1-hour averaged data with the exception of three runs, run 7 on !1 
Oct. and runs 6 and 7 on 12 Oct., which are 40-minutes averaged data. In examining the 
bottom and surface boundary layers, the current meter data are treated separately. For the 
bottom boundary layer, logarithmic profiles based on a least square method were fit to the 
data to test the validity of the logarithmic profile hypothesis (Fig.5). The linear correlation 
coefficients ranged from 0.89 to 0.99 (Table 1). Since the logarithmic velocity profile 
hypothesis is a bottom boundary layer concept, information from current meters near the 
surface can be contaminated by the influences of waves and wind and coming in and out of 
the water and, hence were not included. A criterion is used such that only information from 


current meters below (MSL- H__), where MSL stands for mean sea level, is considered to 
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define logarithmic profiles of the bottom boundary layer. This criterion assures the current 
meters in the analysis came out of the water less than 0.25% of the time based on a eras 
distribution, which is conservative for the measured skewed distributions. The rms wave 
height is approximated by H_ = 80° where o° is the variance calculated from the surface 
elevation time series. The surface elevation was calculated by Fourier transforming the one- 
hour pressure record, applying a linear wave theory transfer function to the complex Fourier 
amplitudes in the frequency domain, and inverse transforming to obtain the surface elevation 
time series. 

In the surface boundary layer, the upper current meters sometimes came in and out 
of water, which causes noise in the current meter output. To eliminate this noise, when the 
current meters were within 5 cm beneath the surface or higher as determined from the 
pressure sensor time sertes, the current velocities were set to zero. 

An example of combining the surface and bottom boundary layers is given in Fig. 
6. The bottom is described by a logarithmic profile and the surface boundary layer by the 
measured surface elevation pdf (equation 11). 

The velocity profiles at successive offshore positions (runs) that the sled occupied 
during a transect are shown in Fig. 7. The data agree well with the logarithmic profile. The 
largest discrepancies occur over the bar, where wave breaking is most intense (x-distance 
between 220 and 240 m). The largest deviation occurs for the fourth run of 11 Oct, which 
generated an outlier in the parameter estimates. The average correlation for all uw profiles 
between log elevation and longshore velocity is 0.98; this value is decreased by the 


observations in the breaking zone (over the bar). 
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The value of z, is calculated from the intersection of the linear regression least 
squares fit to the data in semi-log form, and the shear stress velocity v, 1s estimated from 
the slope. In addition, the bed shear stress coefficient, C pp call be calculated from the data. 


Assuming a quadratic bed shear stress relationship 


1 


ty = pc, (u?+ v)2y (28) 
and combining with (8), gives 
2 
(u 24 y2yi2 v 
A mean eddy viscosity for each run is defined by integrating (10) 
v= bv = (30) 


Values of z,, v,, C, and v, are given in table 1. 

An attempt is made to find relationships between the calculated parameters C,, z, 
and v, with physical measurements such as = (ratio of near-bottom wave velocity 
magnitude to mean current measured at lowest current meter), y= a (Y% of wave 
breaking) and the rms bottom roughness, 7, measured by the CRAB’s altimeter. An example 
of the variation of rms bottom roughness calculated over 20 m detrended sections, in one 
meter steps across the surf zone is shown in Fig. 8. Bottom roughness is smoothest offshore 
and over the bar with increased roughness within the trough. 
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The data are qualitatively sorted by location over the bar, in the trough and on the 
foreshore. This method allows a better identification of the possible soclacone among 
variables. 

In general, C . and z_ values increased an order of magnitude across the surf zone. C f 
values were lowest over the bar and tended to vary over a small range for all values of vm 

y 
(Fig. 9). C , tended to increase with increasing 2 over the trough area. = tended to 
behave in a similar manner to C : (Fig. 10). C ’ does increase with bottom roughness 
(Fig.11). C, is also plotted against — (Fig. 12) where A is the wave orbital excursion 


amplitude, which for linear shallow water wave theory is given by 


A. = PLS y 71 nee (31) 


where g is the gravity constant. = also tends to increase with 7 Fig 13). 

Analyzing the parameters over the cross-shore distance, C , increases slightly toward 
the beach but this smooth trend 1s interrupted with a sharp peak in the wave breaking area. 
After this C_ returns to the previous behavior (Fig. 14). v, shows minimum values over the 


bar, a peak in the trough region and decreasing values on the foreshore, reflecting that it is 


Arms 





Zz 
calculated as the product of v,h(Fig.15). Plots of C 7 oF 7 versus Y= show poor 
correlation (Figs.16 and 17). 
The wind contribution to the longshore currents was evaluated by plotting an 


exponential approximation (26), used to model the wind-driven portion of the profiles, and 


compared with the residual of the measured surface elevation pdf (Fig.18). The exponential 
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approximation shows, in general, good agreement with the residual. The wind data is 


summarized in Table 2. 
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V. DISCUSSION 


There has been a lack of field data to describe the mean vertical velocity profile of 
longshore currents. Three days of DUCK94 experiment when strong longshore currents 
occured were examined. Mean currents obtained using seven current meters spaced from 41 
to 257 cm above the sea bottom, were well described by a logarithmic profile. The validity 
of the logarithmic approach is reflected in a high correlation mean value of 0.98 for all 22 
profiles, between the log elevation of the current meters and the measured velocities. The 
largest deviations occured over the bar and the foreshore, related to the increased turbulence 
caused by wave breaking. 

A sensitivity test was performed to determine if the correlation coefficient changed 
with the number of current meters used to define the profile. Using either all the current 
meters available under the criteria established in chapter 4, or only the lowest three or four, 
the correlations are similar with maximum changes of only 0.02. This demonstrates the 
validity of the logarithmic approach over most of the water column. In the same way, the 
values of C , do not significantly change when calculated with only the lowest current 
meters. 

In the breaking zone the values of the bed shear stress coefficient (or bottom friction 
coefficient), C,, are in the order of 10~* while the rest of the transect shows C, values in the 


order of 10° (consistent with previous experiences in the same beach). 
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There is some uncertainty of the exact distance of the current meters from the bed, 
because the sled runners sunk into the sand an unknown amount depending on bearing 
capacity of the bed (theory and divers suggest 3-8 cm) and the undulating bottom, 
particularly over mega-ripples. To test the sensitivity of C ; due to the uncertainty in depth, 
the elevation of the current meters was shifted +20cm in steps of 1 cm and C , recalculated. 
An average value of aC 7 = 0.001/10cm was calculated (5 to 30 % of calculated C ay The 
changes of C , are less significant for increased depth (-20 cm) than for decreased depth 
(+20cm). 

Z u 

Cand a show surprisingly poor correlation with the wave/current parameter ee 
since it was expected an increasing of C - and with increasing values of -* (increasing 
of the wave influence over the vertical profile). Equation (6) suggests a relationship between 


z, and wave parameters . For shallow water waves (6) can be rewritten 


Z H 
a & ace 0.016 Je rms ae 





(32) 


Equation (32) combines wave, current and roughness variables. 7 is plotted against 
equation (32) expression in Fig.19. However, including wave and current variables with 
roughness makes the correlation withz, worse. 

The strongest correlation was found betweenC ¢ and the bottom roughness, with C - 
increasing with the bottom roughness. Theoretically this is expected as larger roughness 
implies larger bottom stress due to form drag, and consequently larger C , values. Zz, also 
increased with the bottom roughness (-), but no as strongly as C The C . values varied 
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by an order of magnitude across the surf zone. This implies that the bottom roughness has 
an important role in the characterization of the bottom boundary layer by affecting the value 
of C.. Roughness values increase toward the beach and are larger in the trough. The 
increased roughness of the bottom in the trough is associated with mega-ripples observed 
by side-scan sonar (Swayne, 1995). 

The behavior of v, along the cross-shore line is basically a function of the velocities 
and the depth. Low values over the bar and foreshore (smaller / ) and high values over the 
trough (larger velocities) are explained by (30). 

Wind contribution to the longshore current was examined. Using a simplified 
approach, reasonable agreement was found between an exponential approximation and the 
residual of velocity/surface elevation pdf correction. Nevertheless it is difficult to predict if 
the residual is due only to the wind or if other factors, such as mass transport velocity or 
wave breaking also contribute to this residual. The free parameter B in equation (26) varied 
between 3 and 4 — reasonable predictions for the cases studied. Since the approximation 


was tested over only 12 profiles, future work is needed to refine the value of f. 


Zi 
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VI. CONCLUSIONS 


The vertical structure of mean longshore currents on a.barred beach is well described 
by a logarithmic profile for the three strong longshore current days examined. This 
hypothesis works better in the trough where turbulent bottom boundary layer processes are 
more dominant than over the bar, where breaking-wave induced turbulence generated at the 
surface modifies the profile. 

A relationship between near bottom vertical velocity profiles and bottom roughness 
was found. While the relatively small amount of data used (three days) are not enough to 
establish a definitive conclusion from the available data, it can be concluded that the bed 
shear stress coefficient C 71s directly proportional relationship to the bottom roughness and 
therefore is an onan parameter to characterize the bottom boundary layer. Surprisingly, 
only poor correlation of C , and z, was found with wave parameters. 

The influence of the wind is parameterized by an exponential approximation that 
works well in 50% of the analyzed profiles. This parameterization assumes that the residual 
between the velocity corrected using the surface elevation pdf and the measured data is only 
due to the wind and does not consider the alongshore component of mass transport velocity 
that is assumed small. 

These results should be treated carefully since there were obtained from a small 
number of days over a barred beach with strong longshore currents and cannot be 


generalized without further studies. 
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APPENDIX A. FIGURES 
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Figure 1. SLED on the beach during DUCK 94 experiment. 
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Figure 2. Bathymetry and equipment positioning during DUCK 94 experiment (11 Oct.). 


Sled transect line was 940 m alongshore. 
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DUCK 94 - Bathymetry profiles 
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Figure 3. Bathymetry profiles for the three days considered (10, 11 and 12 Oct.). 
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Figure 4. Climatology for the three days considered (10, 11 and 12 Oct.). Currents were 


measured in the middle of the trough. H is the significant wave height and T is the period 
of peak frequency. 
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DUCK 94 - Log Longshore current - Oct 10 
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Figure 5. (a) Least-squares fit to the data for 10 Oct. 
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DUCK 94 - Log Longshore current - Oct 11 
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Figure 5. (b) Least-squares fit to the data for 11 Oct. 
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DUCK 94 - Log Longshore current - Oct 12 
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Figure 5. (c) Least-squares fit to the data for 12 Oct. 
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DUCK 94 - Longshore current - Oct 11 time 15:44 - xposition = 157.3 
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Figure 6. Example of vertical profile (the upper two current meters sometimes were out 
of the water). 
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Figure 7. (a) Vertical profiles of mean longshore currents for 10 Oct. 
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Figure 7. (b) Vertical profiles of mean longshore currents for 11 Oct. 
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Figure 7. (c) Vertical profiles of mean longshore currents for 12 Oct. 


43 











DUCK94 - Bottom roughness - Oct12, Line 10 
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Figure 8. Variation of the Rms Bottom Roughness with the Cross-Shore Distance 
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Figure 9. Relationship Between C_ and U / V. 
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Figure 10. Relationship Between - and U_ / V. 
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DUCK94 - Cf vs Bottom roughness 
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Figure 11. Bed Shear Stress versus Bottom Roughness. The Line Represents a Linear 
Regression Based on All Points but the Outlier. 
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Figure 14. C, versus Cross-Shore Distance. 
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DUCK94 - Eddy Vis. vs X-distance 
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Figure 15 v, versus Cross-Shore Distance. 
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DUCK94 - Cf vs Gamma 
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Figure 16. Relationship BetweenC. and y 
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Figure 17. Relationship Between v, and y 
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DUCK 94 - Longshore current Oct 12 time 11:2 xposition = 251.5 
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Figure 18. Example of wind contribution to the longshore current: exponential 
approximation to the residual. 
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DUCK94 - Za/h vs Hrms r/ h4(3/2) v* 
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Table 1. Profile Fitting Results. 
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Table 2. Wind Data 
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